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Abstract 

We calculate electric and baryonic charge fluctuations on the lattice. Results have been obtained 
with the highly improved staggered quark action (HISQ) and almost physical quark masses on 
lattices with spacial extent of N T = 6,8, 12. Higher order cumulants of the net-charge distri- 
butions are increasingly dominated by a universal scaling behavior, which is arising due to a 
critical point of QCD in the chiral limit. Considering cumulants up to the sixth order, we observe 
that they generically behave as expected from universal scaling laws, which is quite different 
from the cumulants calculated within the hadron resonance gas model. Taking ratios of these 
cumulants, we obtain volume independent results that can be compared to the experimental mea- 
surements. Such a comparison will unambiguously relate the QCD transition temperature that 
has been determined on the lattice with the freeze out temperature of heavy ion collision at LHC 
and RHIC. 



1. Introduction 

Calculating the phase diagram of strongly interacting matter is one of the most important and 
outstanding problem of non-perturbative QCD. A generic plot of the QCD phase diagram based 
on model calculations and model independent symmetry arguments is shown in Fig[T] Due to the 
notorious sign problem lattice QCD calculations are currently restricted to zero baryon chemical 
potential. However, unlike in heavy ion experiments, one can vary the quark mass in lattice QCD 
simulations. As one expects a critical point in the phase diagram at nonzero density, fluctuation 
observables are the most natural choice in order to detect QCD critical behavior. Their divergent 
behavior can give insight in the structure of the phase diagram. Currently various experimental 
programs aim at probing the phase diagram at nonzero baryon number density, such as the low 
energy scan at RHIC, NA61 at the SPS as well as future experiments at FAIR and NICA. In 
experiments one has access to the fluctuations of conserved charges at freeze out. The freeze-out 
curve in the phase diagram can be parametrized as a function of collision energy (TJ. 

2. Cumulants of conserved charges 

On the lattice we study derivatives of the partition function with respect to baryon (B), elec- 
tric charge (Q) and strangeness {S) chemical potentials, which are also kwon as generalized 
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Figure 1: Generic phase diagram of QCD, based on model calculations and model independent symmetry arguments. 
Also indicated are the regions in the phase diagram where we are able to obtain results on fluctuation observables from 
lattice QCD and experiments, respectively. 



susceptibilities and are defined as 



with fix - f*x/T and X = B,Q,S . The lattice studies are performed at fi x — and sufficiently 
close to the thermodynamic limit. The generalized susceptibilities are thus intensive quantities. 
They can also be interpreted as Taylor expansion coefficients of In Z and, furthermore, be related 
to the cumulants of the fluctuations of the net charges (Nx), which are measured in heavy ion 
collisions. E.g., for the diagonal fluctuations one obtains 



fluctuations of net baryon number and electric charge, obtained with an highly improved stag- 
gered quark (HISQ) action on lattices with temporal extent of N T = 6 and 8. In oder to compare 
the lattice QCD results with experimentally measured fluctuations one has to eliminate the un- 
known fireball volume by taking ratios of cumulants. Our approach to fix the strangeness and 
electric charge chemical potentials and to extract the remaining freeze-out parameters, i.e. the 
freeze-out temperature (T-f) and the freeze-out chemical potential from a comparison of 
lattice and experimental results of ratios of cumulants of the electric charge fluctuations has re- 
cently been proposed by us [2 1 and was also presented at this conference J3). At nonzero baryon 
density the method is based on a next to leading order (NLO) Taylor expansion of two different 
ratios. In general we find that the NLO contributions are below 10% in the fig range relevant for 
the RHIC energy scan down to collision energies of yfs > 20 GeV. 



(VT 3 ) -4f CO = (ff + ^ k \nZ(T !MB , M Q,fis))l(d^ B dp^dfi k s ) , 
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Figure 2: Cumulants of the net baryon number fluctuations (left) and electric charge fluctuations (middle) up to the sixth 
order; quadratic cumulant of strangeness, baryon number and electric charge normalized by the HRG results (right). 
Different symbols denote different lattice spacings: N T = 6 (triangles), 8 (circles) and 12 (squares). Solid lines on 
the right panel indicate spline fits to the data. For the electric charge fluctuations additional solid lines indicate HRG 
results with a modified pion mass (see text for explanation). The vertical bar indicates the QCD transition temperature 
as obtained in (7). 



3. 0(4) critical behavior 

In the chiral limit of QCD with two massless quarks the transition from the hadronic world 
to the quark gluon plasma is presumably of second order (see Fig. [TJ. At vanishing chemical 
potential the critical point is expected to be in the universality class of the three dimensional 
0(4) symmetric spin model and is expected to persist even in the presence of a physical strange 
quark mass. It is illuminating to perform a universal scaling analysis that is connected with the 
critical point in the chiral limit. Here to leading order the chemical potential enters only in the 
temperature like scaling field 

t = t 1 ((T - T )/T + K B fi 2 B + K S fi 2 s + K BS {i B {i s y (5) 

as a finite chemical potential does not alter chiral symmetry breaking. It is thus easy to see that 
at fix — the contribution from the singular part of the partition functions to the generalized 
susceptibilities follows the pattern 0), 

Xf 2n) ~ \t\ 2 - n ~ a , (6) 

where a is the critical exponent of the specific heat, which is small and negative. The fourth 
order cumulants will thus develop a cusp in the chiral limit, whereas the sixth order cumulants 
are divergent with amplitudes that have different signs below and above T c , The lattice data is 
qualitatively consistent with that picture as can be seen in Fig. [2] (left and middle). A more de- 
tailed scaling analysis using a recent parametrization of the <9(4)-scaling function of the specific 
heat f5] is work in progress. It will allow to obtain various non-universal normalization constants 
that map QCD to the universal (9(4)-symmetric theory such as To,tQ,ho,Kx,Kxy- Some of them 
are of immediate interest, such as the transition temperature in the chiral limit Tq or the curvature 
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k x that characterizes the change of the transition temperature in the direction of the chemical po- 
tential fi x - Recently the HotQCD Collaboration and the BNL-Bielefeld Collaboration performed 
a scaling analysis of the chiral condensate, the chiral susceptibility ||6l|7l and a particular mixed 
susceptibility [8 |. That also determined the normalization constant zo = hJ^ S /to, where /J and 6 
are critical exponents, which fixes the quark mass dependence of the QCD transition (at fix = 0). 
This was used by HotQCD in oder to determine the crossover temperature T c — 154(9) MeV [7| 
at physical quark masses. 

4. Comparison with the Hadron Resonance Gas 

Below the QCD transition temperature we can compare our lattice results with the statistical 
Hadron Resonance Gas (HRG) model, which describes the hadronization process in heavy ion 
collisions quite successfully Q. For the second order cumulants of baryon number, electric 
charge and strangeness fluctuations this comparison is performed in Fig. [2] (right) ifTUl . Here we 
plot the lattice data normalized by the corresponding HRG results. Light blue bands indicate the 
continuum extrapolations based on lattices with temporal extent N T = 6,8, 12. We find that the 
continuum extrapolations of x\ an d^^ approach the HRG from below and are in agreement with 
the HRG at temperatures up to T ~ (150 - 160) MeV. The strangeness fluctuations, however, 
seem to overshoot the HRG and eventually approach the HRG at much lower temperatures. Note 
that for the electric charge fluctuations additional solid lines at low temperature indicate HRG 
results with a modified pion mass. We have chosen the pion masses such that they agree with the 
averaged pion mass (root-mean-square mass ra^ MS ) of the pion spectrum on the lattice at given 
N T . At finite N T the pion spectrum of staggered fermions on the lattice achieve an unphysical 
splitting which leads to an increased m™. In general, we find that the modified HRG provides 
a good approximation to the electric charge fluctuations below T c . 

A comparison of the HRG with higher order cumulants is at present not very meaningful as 
we do not have a continuum result yet. Moreover, for the electric charge fluctuations, which can 
in principle be immediately compared to the experimental results, the distorted pion spectrum 
on the lattice becomes increasingly problematic, as the higher order cumulants are increasingly 
sensitive to an increased m^ MS . One thus has to use even finer lattice spacings such as A^ T = 16 
in order to control this systematic effect. In fact, we find that the HRG results for^^ with a pion 
mass corresponding to the mj™ of our N T - 6 and 8 lattices are indistinguishable from that of a 
HRG with an infinitely heavy pion mass. 
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